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ABSTRACT 

In  compressible  turbulence  models  it  is  assumed  that  the  Favre-mean  velocities  are  suit¬ 
able  approximations  to  the  Reynolds-mean  velocities  in  order  to  close  unknown  terms.  This 
neglects,  in  the  mean  momentum  and  energy  equations,  the  contribution  to  the  stress  and 
work  terms  by  the  mean  of  the  fluctuating  Favre  velocity,  a  quantity  proportional  to  the 
turbulent  mass  flux.  As  the  stress  and  work  terms  do  not  introduce  any  new  unknown  cor¬ 
relations  requiring  closure  in  either  k  —  e  or  Reynolds  stress  closures  and  because  the  exact 
form  of  the  terms  can,  with  little  additional  work,  be  carried  there  is  no  need  to  make  any 
modeling  assumptions.  In  the  Reynolds  stress  equations  the  viscous  terms  appear  naturally 
in  Reynolds  variables  while  the  problem  is  posed  in  Favre  variables.  In  the  process  of  split¬ 
ting  the  viscous  terms  into  the  viscous  transport  terms,  carried  in  Favre  variables,  and  the 
dissipation  terms,  carried  in  Reynolds  variables,  important  contributions  from  the  mass  flux 
appear.  The  accurate  accounting  of  these  terms  is  important  for 
modeling  and  the  retention  of  the  mass  flux  terms  is  important 
turbulent  flows. 
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Introduction 

The  stress  and  work  terms  in  the  mean  momentum,  mean  energy  and  the  Reynolds 
stress  equations  are  usually  modeled  by  assuming  that  the  Favre  mean  velocity  is  a  suitable 
approximation  to  the  Reynolds  mean  velocity.  This  neglects  the  contribution  of  the  mean 
of  the  fluctuating  Favre  velocity  which,  related  to  the  turbulent  mass  flux,  quantifies  the 
difference  between  the  Reynolds  and  Favre  mean  velocities  of  the  mean  flow.  As  the  stress 
and  work  terms  do  not  introduce  any  new  unknown  quantities  in  k  —  e  or  Reynolds  stress 
models  (an  eddy  viscosity  expression  for  the  mass  flux  is  typically  carried)  and  because 
they  can  be  carried  exactly  with  little  additional  complexity  there  is  no  need  or  justification 
for  any  modeling  assumptions.  The  present  article  derives  the  exact  equations  for  the  mean 
flow  and  the  Reynolds  stresses  for  a  compressible  turbulence  keeping  the  neglected  difference 
between  the  Reynolds  and  the  Favre  mean  velocities.  In  so  doing  it  is  hoped  to  make  two 
crucial  points  clear:  1)  that  the  retention  of  the  mass  flux  terms  in  the  several  places  it 
appears  in  the  mean  momentum,  mean  energy  and  Reynolds  stress  equations  is  essential  to 
the  prediction  of  any  nominally  complex  engineering  flows  and  2)  that  a  careful  distinction 
between  Reynolds-averaged  and  Favre-averaged  variables  must  be  made  to  properly  pose  the 
near-wall  problem  and  insure  that  the  near-wall  asymptotics  are  carried  out  consistently. 

In  this  exposition  upper  case  letters  will  be  used  to  denote  mean  quantities  except  in  the 
case  of  the  mean  density,  <  p  >,  as  p  has  no  convenient  upper  case  form.  The  expectations 
will  be  indicated  using  the  angle  brackets  for  time-means,  eg.  <  u,n7  >,  and  the  curly  brackets 
for  the  density-weighted  or  Favre-means,  eg.  {v.Uj},  where  =<  p*n,u;  ><  p  >_1  and 

the  asterisk  denotes  the  full  field,  P*=<P>  +/>'.  The  dependent  variables  are  decomposed 


according  to 

~  Ui  "f-  Uj 

where  <  u,  >=  0 

U*  =  Vi  +  Vi 

where  {uj}  =0 

P*  =  <  P>  +  p' 

where  <  p'  >=  0 

p*  =  P  +  p 

where  <  p  >=  0. 

As  both  the  Reynolds-mean  and  the  Favre-mean  velocities  appear  in  the  evolution  equations 
for  a  compressible  turbulence  it  is  uecessary  to  carry  both  the  Favre  and  the  Reynolds 
decompositions  of  the  velocity  field.  The  mean  of  the  fluctuating  Favre  velocity,  <  v,  >, 
characterizes  the  difference  between  the  Favre  mean  velocity  and  Reynolds  mean  velocities, 
Vi  and  (/;,  as  well  as  the  difference  between  the  instantaneous  fluctuating  velocities: 

Ui  =  Vi  +  <  Vi  > 

Ui  =  Vi  —  <  Vi  >  . 

The  Favre  fluctuation  mean,  a  first-order  moment,  is  related  to  the  turbulent  mass  flux,  a 
second  order  moment,  by 

<  pvi  >=  —  <  p  ><  v,  > 
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which  follows  from  the  definition  of  the  Favre-average  of  the  Favre  fluctuation,  { v, }  =< 
p*vt  ><  p  >-1=  0.  Thus,  apart  from  a  scaling  by  the  local  mean  density,  the  mean  of  the 
fluctuating  Favre  velocity  and  the  mass  flux  are  equivalent  quantities.  The  primes  on  the 
fluctuating  density  have  been  dropped. 

Mathematically  <  Vi  >  represents  the  difference  between  unweighted  and  density-  weighted 
averages  of  the  velocity  field  and  is  therefore  a  measure  of  the  effects  of  compressibility 
through  variations  in  density.  It  plays  an  important  role  in  parameterizing  the  anisotropic 
effects  of  compressibility  associated  with  the  mean  dilatation  and  the  mean  density  gradi¬ 
ents.  Experimentally  it  is  an  important  and  essential  quantity  that  allows  numerical  results 
computed  in  Favre  averaged  variables  to  be  related  to  experimental  results  computed  in 
Reynolds  variables.  Additional  results  as  well  as  an  equation  for  the  mass  flux  have  been 
given  in  Ristorcelli  (1993).  The  present  investigation  focuses  on  the  role  the  Favre  fluctuation 
mean  plays  in  the  the  mean  momentum,  mean  energy  and  the  Reynolds  stress  equations  and 
shows  that  the  current  practice  of  neglecting  it  is  unnecessary,  inconsistent  and  inadequate. 

The  mean  flow  equations 

Substituting  in  the  Favre  decomposition  into  the  Navier-Stokes  equations  and  time  av¬ 
eraging  produces,  without  approximation, 

(<  P  >  Vi),t  +  (<  P  >  VpVi),P  =  - P,i  +£uu  ( U )  -(<  p>  {u,vp}),p 

where  £,>(f/)  =<  p  >  +Uj,i  —  2/3Uq,q £,;].  Correlations  between  the  fluctuating  viscos¬ 

ity  and  the  fluctuating  velocity  have  been  neglected.  Note  that  the  viscous  terms  are  given 
in  terms  of  the  time-averaged  mean  velocity  while  the  problem  is  solved  in  terms  of  the 
Favre-averaged  mean  velocity.  They  are  related  by  Ui  =  Vi  +  <  Vi  >.  The  usual  assump¬ 
tion  used  to  close  the  viscous  term  is  that  Ui  ~  Vi.  Thus  even  the  first-order  equation  is 
modeled  reducing  the  accuracy  of  the  method  at  the  very  lowest  order  in  the  very  region  of 
the  most  practical  (aerodynamic)  interest.  Part  of  the  appeal  of  k  —  e  or  Reynolds  stress 
turbulence  models  in  the  incompressible  turbulence  is  that  the  equations  for  the  mean  flow 
are  carried  exactly:  this  can  also  be  done  for  the  compressible  turbulence.  Substituting 
Oi  =  \'i  -f  <  Ui  >  produces  the  exact  equation  for  the  mean  momentum  equation: 

(<  p  >  Vi),t  +  (<  p>  VpVi),p  =  ~P,i+  [£,J+  <  <rtj  >j,y  -(<  P>  {t\»>p}),p 

where  =<  p  >  [V^,y  +V},i  -2/3 and  <  atJ  >=<  p  >  [<  ?>,  +  <r  v.  >.,  — 2/3  <" 

vq  b%j\-  The  mean  momentum  equation  is  now  carried  exactly.  Note  that  the  additional 
terms,  <  cr^j  >,  do  not  contribute  to  the  mean  flow  equations  if  the  turbulence  is  either 
homogeneous  or  isotropic.  In  the  near- wall  region  where  the  viscous  terms  are  important 
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and  where  there  are  gradients  in  the  mean  density,  the  <  vt  >  are  important  and  need 
to  be  carried.  Data  from  Ma  =  4.5  DNS  computations  of  Dinavahi  and  Pruett  (1993)  in 
unidirectional  developing  wall  bounded  flow  indicate  that  the  approximation  of  f/t  ~  in 
the  wall  bounded  flow  is  inadequate.  This  is  a  nominally  simple  flow,  in  comparison  to 
those  of  practical  interest,  in  which  there  is  a  four-fold  variation  of  the  mean  density  over 
the  boundary  layer.  In  data  taken  from  that  simulation,  shown  in  Figure  1,  it  was  quite 
unexpectedly  found  that  <  v2  >  is  larger  than  either  U2  and  V2.  It  is  large  enough  to  cause 
V 2  and  V2  to  have  different  signs.  This  indicates  that  the  mean  fluid  particle  transfer  is 
in  a  direction  opposite  to  the  net  momentum  transfer  in  flows  with  mean  density  gradients. 
Moreover,  the  contribution  of  the  Favre  fluctuation  mean  to  the  total  viscous  stress  was  found 
to  be  as  large  as  one  third  the  contribution  of  the  Favre  mean  viscous  stress,  in  the  near 
wall  portions  of  the  turbulent  boundary  layer.  In  flows  with  separation  and  re-attachment 
the  second  derivative  of  the  mean  flow  vanishes  leaving,  in  the  viscous  terms,  the  second 
derivative  of  the  Favre  fluctuation  mean.  Clearly  this  will  be  an  important  term  when  there 
is  a  turbulent  mass  flux  due  to  mean  density  gradients  near  the  point  of  separation. 

The  mean  energy  equation  is  subject  to  similar  deficiencies.  The  exact  equation  for  the 
mean  total  energy,  internal  plus  kinetic,  is 

(<  P  >  E),t  +  (<  P  >  VpE),p=  -  [PVp  +  P  <vp>  +  <pvp  >],p  T  [<  (TpfcUfc  >  -  <  p  >  {eup}],p 

+  [SpfcI4-F  <  o-pjfc  >  14  +  Spit  <  vk  >]>p  QP,p- 

It  is  not  unusual  in  compressible  turbulence  models  to  see  several,  if  not  all,  of  the  terms 
involving  the  mass  flux,  P  <  vk  >,  <  crp*  >  V*,  or  Y,pk  <  vk  >  dropped  because  the  present 
models  for  these  terms  destabilize  computations  or  because  they  are,  in  the  spirit  of  f/,  ~  V,, 
assumed  negligible. 

The  Reynolds  stress  equations 

The  second  moment  equations  for  a  compressible  flow,  are  written  without  approxima¬ 
tion,  and  after  some  manipulation,  as 

D/Dt  (<  p  >  =  -  <  p  >  {viVp}Vjip  -  <  p  >  {vjVp}Vi,p  +  +2/3  <  pvk,k  >  6{j 

-[<  pVi  >  6pj  +  <  PVj  >  btp+  <  p  >  {viVjVp}-  <  Vj(T?p  >  -  <  vta*p  >],p 

+  <  Vj  >  [~P,i+  £;*,*+  <  <rik  >,*]+  <Vi>  \-P,j+  Zjidc  +  <  <y]k  >,fc] 

-  <  Uj,p<7 £  >  -  <  U„p<7",  > 

where  the  mean  momentum  equations  have  been  used  and  cr£  =<  p  >  [uUJ  +Uj»  -2/3 u,„  £,->]. 

The  form  of  the  equations  above  reflect  the  following  manipulations:  1)  The  Jcviatoric  part 
of  the  pressure-strain  correlation  is  defined  f].;  =<  P(voi  +u>>: )  >  —2/3  <  pvk,k  >  8iy  and 
2)  the  identity  vx  —  u;  +  <  Vi  >  has  been  used  to  rewrite  the  transport  terms  in  Vi  variables 
while  keeping  the  dissipation  terms  in  U{  variables.  In  the  Reynolds  stress  equations  the 
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terms  arising  from  viscous  surface  terms  appear  naturally  in  u,  variables  while  the  problem 
is  posed  in  the  Favre  u,  variables. 

The  dissipation  type  terms,  which  require  closure,  are  kept  in  ux  variables.  The  fact 
that  Vi  =  «,■  +  <  V{  >  is  used  to  segregate  the  terms  into  ones  that  require  modeling 
and  ones  that  are  carried  as  dependent  variables  in  the  closure  scheme.  The  identities 

<  m,pup,i  >  =  <  UiUp  >,ip—  2  <  Ui ,,  up  >,p+  <  Ui,iUp,p>  and  Ui,p  =  tip„  -  eqpco,  wh^re  the 
vorticity  is  u>i  =  Eijk^jtk  are  used  to  write  the  trace  of  the  dissipation  terms  as 

^  Uj  »P  ^ jp  ^  ft  ^  ^k^k  "t"  2(  XI j Up  'P*  xjp  2  Ukxk  Up  ^ ,p  )  +  4/3  Xlkxk  Up,p  >]. 

Here  u)jU)j  ^ Uj,ptij,p^>  UjxpUpxj  UjxpUjxp^>  UgUp  ^>xqp-\~2  uqxqup 

,p  —  <  uq,gUp, p>.  In  a  homogeneous  compressible  turbulence  <  Wju;;  >  =  <  u],Pu],v>  —  < 
Uj,jUp,p>  which  reduces  to  the  usual  expression  in  an  incompressible  homogeneous  turbu¬ 
lence.  Defining  the  positive  definite  scalar  dissipation  quantities,  the  solenoidal  dissipation 

<  p  >  e3  =<  p  ><  u ijCJj  >  and  the  dilatational  dissipation  <  p  >  Ed  =  4/3  <  p  >< 
Up,puq,q  >  the  trace  can  then  be  written  as 

<  Uj,p  (Tjp  >=  <  p  >  (e,  +  £d)  +  et 

where  et  =  2  <  /i  >  [<  UjUp  >,7  —2  <  u*,*  up  >],p  is  a  scalar  transport  term  that  comes  from 
the  dissipation  type  terms.  The  decrease  of  the  Favre  kinetic  energy,  for  a  homogeneous 
turbulence  is  then  written  as  D/Dt  (<  p  >  {ujVj})  =  —2  <  uJxpcru3p  >—  —2  <  p  > 
(e*  +  £d)  — 2et  where  et  =  0  Note  that  st  is  not  defined  as  per  unit  mass  quantity  to  emphasize 
that  an  equation  need  not  be  carried  for  it:  the  substitution  u,  =  Vi  —  <  Vi  >  shows  that 
it  can  be  written,  except  for  the  correlation  with  the  divergence,  in  terms  of  the  mass  flux 
and  the  Reynolds  stresses  for  which  equations  are  carried.  In  general  et  is  either  positive 
or  negative:  in  a  homogeneous  turbulence  et  is  zero  while  in  the  near-wall  region  it  makes 
a  nonnegligible  contribution  to  the  energy  budget.  The  point  is  that  in  an  inhomogeneous 
turbulence  it  is  necessary  to  recognize  the  contributions  to  <  >  that  are  functions  of 

gradients  of  the  Reynolds  stresses  rather  than  hiding  these  terms  in  the  dissipation.  If  the 
trace  is  added  and  subtracted  the  dissipation-type  terms  can  be  rewritten  as 

<  Wj,p(Tjp  >  <  «,\p  crjp  <  p  >  2/3  <  p  >  (e 3  T  £d)8ij  2/3EtSij 

where  the  term  <  P  >  4  has  zero  trace.  The  viscous  terms  are  now  manipulated  into  their 
final  form.  Using  cr,"  =  <rq  —  <  cq,  >,  which  follows  from  u;  =  Vi  —  <  v,  >,  allows  the 
viscous  transport  terms  to  be  rewritten  as  <  t>icr“p  >  =  <  v,a]p  >  —  <  u,  ><  erJP  >.  The 
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Reynolds  stress  equations  can  be  written  in  the  following  form 
D/Dt  (<  p  >  {v^j})  =  -  <  p  >  {viVp}Vj,p  -  <  p  >  {vjVp}V„p  +  Y\tJ  +2/3  <  pd  >  StJ 

-[<  pvt  >  6pj+  <  pvj  >  8ip+  <  p  >  {vtVjVp}],p 

+  2[<  p  ><  v}Vi  >,p  -  <  p  ><  Vid  >  8jp/ 3-  <  p  ><  Vjd  >  6xp/ 3],p 

-2/3 et8ij  -  [<  Vj  ><  crip  >  +  <  Vi  ><  aJP  >],p 

+  <C  Vj  >  [  Pyi  +  E,p,p  +  <  (7 ,p  >?p]+  ^  +  Yjjp,p  +  <  <7jp  >,p] 

-  <  p  >  ej  -  2/3  <  p  >  (e,  +  ed)8tJ 

where  the  definition  of  the  stress  tensor  <r,j  has  been  used  to  expand  [<  Uj<T,p  >  +  <  Vi<Jjp  >] 
in  the  viscous  transport  terms  and  d  =  vp,p.  Except  for  neglecting  correlations  between  the 
fluctuating  Favre  velocity  and  its  vorticity  the  final  form  of  the  Reynolds  stress  equations, 
given  above,  is  exact. 

Conclusions 

The  mass  fluxes  or  Favre  fluctuation  means  appears  in  several  places  in  the  compressible 
turbulence  equations.  They  contribute  to  1)  the  viscous  terms  in  the  mean  momentum 
equation,  and  in  the  mean  energy  equation  they  contribute  to  2)  the  viscous  terms,  3)  the 
pressure  work  terms,  4)  the  viscous  work  terms,  and  the  5)  pressure  flux  which  is  coupled  to 
the  mass  flux  through  the  equation  of  state.  Modeling  f/,  ~  K  is  equivalent  to  neglecting  the 
mass  flux  in  five  different  locations  in  the  mean  equations.  In  the  Reynolds  stress  equations 
the  mass  flux  contributes  to  6)  the  viscous  diffusion  of  the  Reynolds  stresses,  which  only 
recognized  when  the  viscous  terms  are  properly  distinguished  into  their  Reynolds  and  Favre 
variable  components.  The  mass  flux  determine  the  importance  of  two  production  mechanisms 
7)  one  due  the  acceleration  of  the  mean  flow  and  8)  the  other  due  to  viscous  effects  associated 
with  the  Favre  fluctuation  mean  and  also  it  contributes  to  9)  the  pressure  flux.  A  general 
algebraic  model  for  the  <  >,  derived  from  first  principles,  and  suitable  for  complex  flows 

of  engineering  interest  has  been  derived  in  Ristorcelli  (1993)  and  tested  in  Ristorcelli  and 
Dinavahi  (1993). 

The  viscous  terms  in  the  Favre- averaged  Reynolds  stress  equations  have  been  system¬ 
atically  derived,  identified  and  segregated.  The  source  of  terms  associated  with  anisotropy 
and  inhomogeneity  of  the  dissipation,  essential  to  consistent  near-wall  modeling,  have  been 
shown.  The  second  essential  point  is  that  the  dissipation  is  a  function  of  u,  fluctuations 
about  the  Reynolds  mean,  while  the  Reynolds  stress  equation  is  an  evolution  equation  for 
statistics  of  u,  fluctuations  around  the  Favre  mean.  The  near-wall  Taylor  series  expansions 
of  these  two  quantities,  in  a  general  compressible  flow,  are  different.  Near-wall  asymptotics 
must  recognize  these  facts  when  relating  Favre  type  variables  in  which  the  prob¬ 

lem  is  posed,  to  dissipation  type  quantities  which  are  carried  in  Reynolds  variables.  This 
point  is  crucial  for  any  consistent  general  near-wall  model  development  in  flows  only  nomi- 
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Daily  more  complex  than  the  usual  uni-direction  adiabatic/isothermal  flat  plate  flows:  it  is 
necessary  to  start  with  the  exact  equations  distinguishing  between  Favre  and  Reynolds  de¬ 
composed  variables  and  recognizing  that  their  respective  Taylor  series  expansions  are  related 
by  vj  =  Uj+  <  Vj  >. 

It  is  typical  of  current  compressible  k—e  and  Reynolds  stress  models  to  neglect  some  or  all 
of  these  contributions  to  the  mean  momentum,  mean  energy  and  Reynolds  stress  equations. 
The  retention  of  the  mass  flux  terms  will  be  important  in  complex  compressible  turbulent 
flows:  these  include  flows  in  which  there  are  mean  density  gradients  due  to  large  Mach 
number  or  combustion,  separation  or  reattachment  (inflection  points),  cold  wall  boundary 
conditions,  mean  dilatation,  shocks,  adverse  pressure  gradients,  or  strong  streamwise  accel¬ 
erations. 
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Figure  1.  The  Reynolds  mean,  Favre  mean  and  Favre  fluctuation 
mean  velocity  profiles. 
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